Bayesian Prediction of Severe Outcomes in the LabMarCS: Laboratory Markers of COVID-19 Severity - Bristol Cohort (preprint)

Objectives: To develop cross-validated prediction models for severe outcomes in COVID-19 using blood biomarker and demographic data; Demonstrate best practices for clinical data curation and statistical modelling decisions, with an emphasis on Bayesian methods. Design: Retrospective observational cohort study. Setting: Multicentre across National Health Service (NHS) trusts in Southwest region, England, UK. Participants: Hospitalised adult patients with a positive SARS-CoV 2 by PCR during the first wave (March - October 2020). 843 COVID-19 patients (mean age 71, 45% female, 32% died or needed ICU stay) split into training (n=590) and validation groups (n=253) along with observations on demographics, coinfections, and 30 laboratory blood biomarkers. Primary outcome measures: ICU admission or death within 28-days of admission to hospital for COVID-19 or a positive PCR result if already admitted. Results: Predictive regression models were fit to predict primary outcomes using demographic data and initial results from biomarker tests collected within 3 days of admission or testing positive if already admitted. Using all variables, a standard logistic regression yielded an internal validation median AUC of 0.7 (95% Interval [0.64,0.81]), and an external validation AUC of 0.67 [0.61, 0.71], a Bayesian logistic regression using a horseshoe prior yielded an internal validation median AUC of 0.78 [0.71, 0.85], and an external validation median AUC of 0.70 [0.68, 0.71]. Variable selection performed using Bayesian predictive projection determined a four variable model using Age, Urea, Prothrombin time and Neutrophil-Lymphocyte ratio, with a median AUC of 0.74 [0.67, 0.82], and external validation AUC of 0.70 [0.69, 0.71]. Conclusions: Our study reiterates the predictive value of previously identified biomarkers for COVID-19 severity assessment. Given the small data set, the full and reduced models have decent performance, but would require improved external validation for clinical application. The study highlights a variety of challenges present in complex medical data sets while maintaining best statistical practices with an emphasis on showcasing recent Bayesian methods.

Association of COVID-19 with arterial and venous vascular diseases: a population-wide cohort study of 48 million adults in England and Wales (preprint)

Importance: The long-term effects of COVID-19 on the incidence of vascular diseases are unclear. Objective: To quantify the association between time since diagnosis of COVID-19 and vascular disease, overall and by age, sex, ethnicity, and pre-existing disease. Design: Cohort study based on population-wide linked electronic health records, with follow up from January 1st to December 7th 2020. Setting and participants: Adults registered with an NHS general practice in England or Wales and alive on January 1st 2020. Exposures: Time since diagnosis of COVID-19 (categorised as 0-6 days, 1-2 weeks, 3-4, 5-8, 9-12, 13-26 and 27-49 weeks since diagnosis), with and without hospitalisation within 28 days of diagnosis. Main outcomes and measures: Primary outcomes were arterial thromboses (mainly acute myocardial infarction and ischaemic stroke) and venous thromboembolic events (VTE, mainly pulmonary embolism and lower limb deep vein thrombosis). We also studied other vascular events (transient ischaemic attack, haemorrhagic stroke, heart failure and angina). Hazard ratios were adjusted for demographic characteristics, previous disease diagnoses, comorbidities and medications. Results: Among 48 million adults, 130,930 were and 1,315,471 were not hospitalised within 28 days of COVID-19. In England, there were 259,742 first arterial thromboses and 60,066 first VTE during 41.6 million person-years follow-up. Adjusted hazard ratios (aHRs) for first arterial thrombosis compared with no COVID-19 declined rapidly from 21.7 (95% CI 21.0-22.4) to 3.87 (3.58-4.19) in weeks 1 and 2 after COVID-19, 2.80 (2.61-3.01) during weeks 3-4 then to 1.34 (1.21-1.48) during weeks 27-49. aHRs for first VTE declined from 33.2 (31.3-35.2) and 8.52 (7.59-9.58) in weeks 1 and 2 to 7.95 (7.28-8.68) and 4.26 (3.86-4.69) during weeks 3-4 and 5-8, then 2.20 (1.99-2.44) and 1.80 (1.50-2.17) during weeks 13-26 and 27-49 respectively. aHRs were higher, for longer after diagnosis, after hospitalised than non-hospitalised COVID-19. aHRs were also higher among people of Black and Asian than White ethnicity and among people without than with a previous event. Across the whole population estimated increases in risk of arterial thromboses and VTEs were 2.5% and 0.6% respectively 49 weeks after COVID-19, corresponding to 7,197 and 3,517 additional events respectively after 1.4 million COVID-19 diagnoses. Conclusions and Relevance: High rates of vascular disease early after COVID-19 diagnosis decline more rapidly for arterial thromboses than VTEs but rates remain elevated up to 49 weeks after COVID-19. These results support continued policies to avoid COVID-19 infection with effective COVID-19 vaccines and use of secondary preventive agents in high-risk patients.

Understanding COVID-19 trajectories from a nationwide linked electronic health record cohort of 56 million people: phenotypes, severity, waves & vaccination (preprint)

Background: Updatable understanding of the onset and progression of individuals COVID-19 trajectories underpins pandemic mitigation efforts. In order to identify and characterize individual trajectories, we defined and validated ten COVID-19 phenotypes from linked electronic health records (EHR) on a nationwide scale using an extensible framework. Methods: Cohort study of 56.6 million people in England alive on 23/01/2020, followed until 31/05/2021, using eight linked national datasets spanning COVID-19 testing, vaccination, primary & secondary care and death registrations data. We defined ten COVID-19 phenotypes reflecting clinically relevant stages of disease severity using a combination of international clinical terminologies (e.g. SNOMED-CT, ICD-10) and bespoke data fields; positive test, primary care diagnosis, hospitalisation, critical care (four phenotypes), and death (three phenotypes). Using these phenotypes, we constructed patient trajectories illustrating the transition frequency and duration between phenotypes. Analyses were stratified by pandemic waves and vaccination status. Findings: We identified 3,469,528 infected individuals (6.1%) with 8,825,738 recorded COVID-19 phenotypes. Of these, 364,260 (11%) were hospitalised and 140,908 (4%) died. Of those hospitalised, 38,072 (10%) were admitted to intensive care (ICU), 54,026 (15%) received non-invasive ventilation and 21,404 (6%) invasive ventilation. Amongst hospitalised patients, first wave mortality (30%) was higher than the second (23%) in non-ICU settings, but remained unchanged for ICU patients. The highest mortality was for patients receiving critical care outside of ICU in wave 1 (51%). 13,083 (9%) COVID-19 related deaths occurred without diagnoses on the death certificate, but within 30 days of a positive test while 10,403 (7%) of cases were identified from mortality data alone with no prior phenotypes recorded. We observed longer patient trajectories in the second pandemic wave compared to the first. Interpretation: Our analyses illustrate the wide spectrum of severity that COVID-19 displays and significant differences in incidence, survival and pathways across pandemic waves. We provide an adaptable framework to answer questions of clinical and policy relevance; new variant impact, booster dose efficacy and a way of maximising existing data to understand individuals progression through disease states.